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The main objective of this 3tudy has been to 
determine whether the principles of linear wave super- 
position could be combined with model tests using wave 
survey techniques to obtain a practical method of 
deriving "influence diagrams’ 1 f which relate the influ- 
ence of small hull form variations on wave resistance. 
Tests were run in the ..'ebb towing tank using a 5-5 ft 
Destroyer-type form with protuberances added at a number 
of locations over a range of speeds. Results showed 
attempts to apply linear superposition for this purpose 
to be largely unsuccessful. This result implies serious 
limitations of the linearized theory of wave resistance 
and its application to develop optimum hull forms. 
However, several of the results point to the conclusion 
that, with some modification, linear superposi tion may 
be successfully applied toward the development of 
influence diagrams. 
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Three wise men of Gotham 

Went to sea in a bowl 

If the bowl had been stronger. 

My story would have been longer. 
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INTRODUCTION 



William w roude*s Investigation of ship resistance of 
nearly one hundred years ago Inferred that waves produced 
by a ship moving through the water are a measure of the 
energy imparted to the fluid. It is only in the past ten 
years, however, that this question has received the close 
attention of hydrodynami cists from an experimental point 
of view. Much of this attention has been focused on the 
theoretical aspects of the wave resistance as derived from 
the wave pattern of a ship, and experimental verification 
by model tests. Some of the most promising experimental 
techniques developed to date have been assessed by Eggers, 
Sharma, and Ward (1)^ with very favorable results. There- 
fore, it is the authors* contention that the state of the 
art has now been advanced to the point where practical 
applications of wave resistance survey techniques may 
be made. 

An Important assumption made in linearized wave theory 
is the purported ability to superimpose the individual wave 
patterns of independent hull form variations to obtain a 
compound wave pattern equivalent to the pattern which would 
have been obtained had the hull and form variations been 
towed simultaneously. 



1. Numbers in brackets indicate references at the end of 
this paper. 
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This paper will report the results of an exhaustive 
Investigation into this linear superposition theory, and of 
attempts to develop it into a practical means of predicting 
the wave resistance trends of Destroyer-type hulls. To 
this end the following technique is Introduced: 

A Destroyer-type model will he tested to determine its 
hare hull wave slope record hy means of the longitudinal 
cut, wave survey method. A protuberance will then be 
attached to the model and the tests replicated. ?rom these 
records the wave resistance of the model may he determined. 
Also, assuming linear superposition of the wave records, 
the wave record of the protuberance alone may he determined 



That is, the wave record of the protuberance will be assumed 
to be unaffected by its location on the model and to give 
the same wave slope record when attached to the model as 
it would if it were alone in a free stream. This hypothesis 
is investigated in depth in a later section of this paper. 

The protuberance wave record will be shifted along the 
model length and added back into the bare hull record. The 
resultant wave record obtained is, by the linear superposi- 
tion assumption, equivalent to the wave record that would 
have been obtained had the model been towed with the 
protuberance at that shifted location on the model. An 
influence diagram, as suggested by Hogner (2) , can be 
developed by considering the protuberance as an increase of 



y protuberance 




model with 

protuberance bare hull 
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the seotional area at the point of attachment. An 
Influence diagram indicates the change in wave resistance 
per unit change of displacement at the different locations 
on the hull surface, and therefore shows at what points 
a given form should be modified in order to reduoe its 
wave resistance. 

Influence diagrams developed completely from experi- 
mental testing require an exorbitant amount of data re- 
cording and analysis. However, if linear superposition 
methods are successful in this paper, diagram construction 
effort will be greatly reduced, and an extremely useful 
tool will be made readily available to the designer. 
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THEORY 



Introduction 

The particular wave measuring technique to be employed 
In this paper will be the longitudinal cut, wave slope 
method (3). which is selected for the following reasons: 

(a) a longitudinal cut through a wave pattern is readily 
obtained by towing a model past a stationary probe, (b) the 
wave slope record terminates much more rapidly than do 
height records, (c) it Is a method recommended by Eggers, 
bharma, and Ward (1,3) » and (d) all of the necessary 
equipment was available at Webb Institute of Naval 
Architecture. 

Basic Equations and Definitions 

The wave theory which forms the basis for the lateral 
wave slope method of survey is developed by Sharma (4) via 
energy flow considerations across a boundary or continuum. 
The final expression relating the wave resistance to its 
respective wave pattern is given by: 




where S v and C are Fourier sine and cosine 

y y 



transforms 



of the free wave slope profile. That is. 




( 2 ) 



( 3 ) 



where 



Is the wave slope corresponding to any given 



value of x, u and w are the longitudinal and transverse 
circular wave numbers respectively, and k Q is a velocity 



required in arriving at a usable solution of equation (1), 
simple wave theory will be reviewed. 



A representative wave of wave length \ (at point A) 
in the wave pattern of a ship advancing with velocity V m 
is sketched In Figure 1. From the axes' orientations it 
can be seen that & is the angle between^) , the resultant 
coordinate of the wave, and the ship velocity. Two 
conditions must be satisfied: namely, that the wave system 
remain stationary with respect to the ship, and that the 
ship be in deep water. Therefore, wave speed is equal to 



letting the subscript zero denote conditions at the ship 
centerline, we have: 



parameter equal 




To fully understand the steps 





With k termed the wave number - 



- and 






c - Vi, COS & -- 4WJ w 

G - V^\fp^[2K -\JfjT 0 (5) 

and A = fjy^ cos*& - J: 0 S£c x & <«> 

Prom Figure 1 It Is seen that (?) 

cj-X oose+Ys/we 



Multiplying equation (6) by equation (7), 

Aa ) =-(? 0 (xcoses£c*e •■ys/we Sec 2 ©) 

A 6) - Aot sue & 4.4y S £ce -fc/w & 

Jico = Y/X +LIV 

Kher.V=-4s£ce f'H=4S£CS*W® 

squaring and adding equations (9)» 

1/ -Vy/V/ 4 -! 



(?) 

(9) 

(10) 
(ID 



The circular wave numbers u and w have thus been expressed 
in terms of known quantities and each other. 

A solution of equation (1) is outlined by Ward (2) 
as follows. The wave numbers can be non-dlmensionalized 

by dividing each by k Q . ^ v'= w /.i 0 = sec© 

A -^ / Ji 0 - secs-mkis 



Equation (1) then reduces to 

•go 



RwirCfv-v) 
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Since k Q = 



( 12 ) 



3/vJ > g W) 

2. The primes will be omitted throughout the remainder of the 



text. 



S Is determined in a straightforward manner where x 

t/ 

varies from zero at some point ahead of the wave to a point, 
N, where the ^ die out. 






( 13 ) 



where each^ corresponds to the wave recorded slope tape 
record. This integration must be carried out for each 
value of u. A plot of (!•«•• the wave amplitude) 
versus u, w, or © gives a measure of the generated energy 
of the wave system. Similar wave patterns will display 
similar wave amplitude spectra. It can be anticipated 
that two wave records possessing seemingly large discrepancies 
may have very similar wave amplitude spectra, because of 
the smoothing effect of the Integration process. A further 
discussion of these spectra will be conducted in a later 
section. 

Equation (12) can be non-dimensionallzed to compute 
the wave resistance coefficient, Cw, since 

R u)~ l /z CwfV^S • C w =f5 y*£w-i) (14) 

The basic equation is now in a final form convenient for 
adaptation of a digital computer program solution. The 
programming details and assumptions will be described in 
Appendix B. 
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EXPERIMENTS 



Characteristics of Model and Protuberances 
Stevens model number 393 was modified for use In this 
project by removing the bilge keels and adding turbulence 
stimulators. Figure 2 Is the model body plan and Table I 
lists the model and ship particulars. The two protuber- 
ances used are shown In Figure 3. Both are bodies of revo- 
lution and, consequently, all sections are circular. 



TABLE I SHIP - MODEL CHARACTERISTICS 



SHIP 

Length of load waterline 
Length between perpendiculars 
Beam 
Draft 

Displacement - tons salt water 
Wetted surface 

c b 

Displacement-Length ratio 
Beam-Draft ratio 



342.5 ft 
341.0 ft 

35.9 ft 
12.7 ft 

2 , 296.0 
14,160.0 ft 2 
0.517 

57.9 
2.33 
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TABLE I 



SHIP-MODEL CHARACTERISTICS (con.) 



MODEL 



Length load waterline 

Length between perpend loulars 

Displacement-pounds fresh water 



Lambda ratio 



62.0 

5.524 ft 



Wetted surface 



5.500 ft 
21.04 
3.682 ft 2 



Measured freeboard aft focsle break 



2.17 In 



Trim 



none 



Number of stations 



10 



Exploratory Testing 

In order to proceed with the proposed objectives, it 
was neoessary to determine approximate towing speeds as well 
as protuberance sizes which would generate significant wave 
slope difference signals. Practical limitations confined 
the size of the protuberances between the design water line 
and the keel. Moreover, selection of a suitable hull location 
was necessary for comparison of the difference signals and 
the shifting technique which would be developed. Regarding 
towing speed. It was felt that the selection should depict 
the typical operational speeds of a Destroyer, since design 
techniques (influence diagrams) were of concern. A three 
wire resistance gage and Brush recording device were used to 
obtain the wave slope records. Equipment details are described 
In Appendix A. 
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Preliminary model runs were made with two hemlspherloally 
shaped protuberanoes attached at station 1 of a 10-station 
model. It was assumed that this forward location would give 
the best representation of a bulb generated wave since it had 
the least hull separation of the protuberanoes . The results 
of these equipment-model familiarization tests tentatively 
verified that the speeds and the selected protuberance size 
would give sufficient difference signals necessary for 
meaningful computer analysis. Consequently, qualitative testing 
could proceed. 

Resistance Tests 

A standard model resistance test was made to obtain 

sufficient data to develop model residuary and total resistance 

curves. The computed values of C^. could therefore be compared 

with C r and C t to determine the relative magnitude of the wave 

resistance at various speeds. It could also serve as a check 

to insure that the computed wave resistance constitutes a 

43 

reasonable percentage of residuary resistance. 

Prom the results of the aforementioned investigations, 
it was decided to tow off-center at Froude numbers of .242, 

. 363 , and .484 (corresponding to ship speeds of 15. 22.5, 
and 30 knots respectively) , with the larger of the two pro- 
tuberances attached successively at stations 1 , 3. and 5. and 



4.^ The authors know of no previous wave resistance testing 
of a Destroyer-type model. Therefore, the expected range 
of values is uncertain. 
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the smaller attached only at station 3. To provide addi- 
tional data for translation analysis techniques, later 
tests were conducted with the large protuberance at stations 
2 and 3» The protuberances would be Investigated at a 
common waterline for comparative purposes. figure 4 shows 
the model being towed at Proud e numbers of 0.242 and 0.484 
with the large protuberance attached at station 2. 

Any errors Introduced when comparing the records at 
each configuration, and the sensitivity of the testing 
equipment must be minimized throughout testing. Additionally, 
the influence on wave resistance of the transverse distance 
between the gage and the model must be consistent. A 
detailed error analysis will be presented in Appendix D. 
Suffice It to say, the following procedures were employed 
to reduce areas of possible error: 

1. A calibration factor relating the recorded values 
in mm to a slope in radians is essential for each tape 
record. Since the equipment gain settings tend to drift 
during the testing, calibration .records were obtained 
prior to and immediately after each testing configuration. 

The gage was physically rotated through known angles and 
the ensuing deflections were recorded on Brush recorder 
tape. The values were plotted to determine an average 
calibration factor for each specific configuration. It 
was noted that these values were a definite function of 
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the equipment settings and that they greatly affect the 
final wave resistance coefficients. 

2. The employment of micro-switches established a 
common wave slope reference point and a rough verification 
of the model speed. One micro-switch was positioned on 
the tow rail at the point where the bow at design water 
line passed the wave slope probe, the other was positioned 
at a point near the end of the run. The switches produced 
a set of marks on each Brush recorder tape record, and by 
aligning these visually, the records at a given Proude 
number could be qulokly compared for consistency. 

3. Bare hull model tests were conducted at the 
middle speed (Fr * . 363 ) with the distance between the 
model and wave probe being varied between 18 and 24 lnohes. 
In order to determine a range of off-set distance where 

C w remained essentially constant and near the maximum. 

The midpoint of this range was utilized throughout all 
towing configurations. Detailed results are shown In the 
data analysis section. 

Three wave slope records were obtained at each Froude 
number for each of seven model-protuberance configurations, 
each speed varying no more than t .01 ft/seo over a 35*foot 
measured run. Figure 5 shows a sample wave record of the 
bare hull model at Proude number .363 and Figure 6 shows 
a sample wave record with the large protuberance located 
at station 3 at Froude number . 363 . It Is worthy of note 
that all testing (total of 63 acceptable records) required 
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FIGURE 6 

SAMPLE WAVE SLOPE RECORD 



PROTUBERANCE AT STATION NO. I 

PR. NO. 0.303 
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only an eleven-hour period. This time Included preparing 
the model for each configuration, reballasting, trimming, 
and calibration of the wave slope measuring probe, as well 
as conducting the series of off-center runs. 
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METHOD OF DATA ANALYSIS 



Since there have been no known previous wave resist- 
ance tests on Destroyer-type models, a method of reading 
and processing the data had to be developed before any 
detailed analyses, comparisons, or protuberance shifting 
could be initiated. Three representative sets of data were 
used to standardize the data reading and processing pro- 
cedures. Data used were for the bare hull configuration 
with the large bulb attached at station 1 at a Froude 
number of 0.242. The following discussion pertains to 
these three sets of data, unless otherwise specified. 

The records obtained from the model tests had to be 
magnified before they could be accurately read. Several 
methods were tried. It was found that the easiest and 
most reliable method was to project the wave slope records 
on the wall next to the computer console and to read the 
data directly onto computer tape. The accuracy of the 
readings thus obtained is estimated to be t 0.2 mm. 

Two alternatives were available for making use of 
the three data tapes obtained for each speed and protuberance 
location. Data tapes could be averaged and a wave resist- 
ance coefficient calculated from the resulting data, or 
a resistance coefficient could be obtained for each data 
tape, and the coefficients averaged for a final wave 
resistance coefficient. Both methods were applied and 
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the difference between wave resistance coefficients 
obtained was stiff lclently small as to be attributable to 
the rounding-off error introduced in averaging the indi- 
vidual data points. However, there was sufficient difference 
between the wave resistance coefficients obtained from the 
individual tapes to preclude use of a single tape as 
representative of the given configuration. 

It was decided to analyze the wave records on the 
basis of an averaged data tape since this method adapts 
itself more easily to shifting the protuberance and requires 
less computer core time. 

Two reading intervals were considered, 2.5 mm and 
5.0 mm. Data tapes for each configuration were read at 
2.5 mm. Intervals. In each case the tape length read was 
600 mm. Wave resistance coefficients were then obtained 
by using (a) all data points obtained (240 points, 2.5 mm. 
interval), and (b) every other point (120 points, 5.0 mm. 
Interval) . Comparison of the wave resistance coefficients 
for a given configuration and speed showed less than 3 a 
deviation, as shown in Table II. It was therefore decided 
to use the 5.0 mm. Interval in further analyses to simplify 
the tape reading procedures. This will have no effect on 
the results since they are comparative in nature. 
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TABLE II 



Have 'Resistance Coefficients at 2.5 and $.0 ram, Reading; Intervals 



Fr. 


Pro tub. 


Fro tub. 


Length 


No. of 


C w at 
2.5 mm. 


C'w at 

5.0 mm . 


% Dlff 




Size 


Loca. 


of Data 


Inter- 


from 








Tape 


vals 


Inter- 


Inter- 


5.0 mm 








Head 


of u 


vals 


vals 


value 






(Sta.) 


(mm. ) 




(xlO 3 ) 


(xlO 3 ) 




0.242 


Large 


1 


240 


100 


0.317 


0.323 


- 1.9 


0.363 Bare 




240 


ICO 


0.648 


0. 666 


- 2.7 




Hull 














0.434 


Bare 

Hull 




240 


100 


2.359 


2.344 


+ 0.6 



Theory states that the integration of wave slope data 
with respect to the transverse circular wave number cover 
the range -<*> to +oo . The basic Interval between wave 
numbers used was 0.1 (3ee Appendix C) and resistance coeffi- 
cients were obtained for lengths of 0 through 125 Intervals 
of u. The results showed that beyond 100 intervals the increase 
obtained in the resistance coefficient was less than 

To determine the maximum number of data points re- 
quired to obtain reliable wave resistance coefficients, two 
procedures were used. First, the data tapes were analyzed 
for lengths varying between 400 and 600 mm. The results 
did not appear to be approaching a limiting value for any 
of the three Froude numbers. Therefore, the tank geometry. 
Figure 7 i was developed, assuming the bow wave angle to be 
equal to the Kelvin angle (19§°)» and a maximum tape length 
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WAVE 



to avoid reflection was determined. The data tapes were 
then re-analyzed. Figure 8 is a graphical representation 
of the final analysis. The final data tape reading lengths 
selected were 800 mm. for a Froude number of 0.242, since 
the wave records die out beyond this tape length, and 700 mm. 
and 550 mm. for Froude numbers of 0.3^3 and 0,484 respectively, 
since reflection occurs beyond these lengths. 

The data obtained for model tests with varying offset 
was analyzed for a data tape length of 600 mm. and a reading 
interval of 5.0 mm. The results, shown as Figure 8, are 
considered to verify that the offset used in the main 
resistance tests (22| inches) was well within the range 
which gives little or no change of wave resistance coeffi- 
cient with change in offset. 
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DISCUSSION OF RESULTS 



The results were compared In a number of different 
ways. These were: comparison of residuary and total re- 
sistance to wave resistance, linear shifts of the difference 
signals, translation of the Influence diagrams, small 
Interval shifts of the difference signals, wave spectra 
analysis, comparison of difference signals, and a compari- 
son of the effects of protuberance size. Each of these 
and representative curves will be discussed In the 
following section. 

Comparison of Residuary and Total Resistance with Wave 
Resistance 

In order to obtain a quantitative comparison of the 
wave resistance component with the total and residual 
components, and thus evaluate the predicted savings In 
resistance In later discussions. Figure 9 has been prepared. 
Table III Is also presented to give a percentage comparison 
between the various resistances. These ratios are for the 
most part In keeping with the results of wave resistance 
tests on merchant type hull forms and are thus considered 
as an indication of the validity of the tests performed 
for this report. 
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BARE HULl MODEl RESISTANCE 
COEEnCiENTS 
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TABLE III 



Resistance Coefficients 


5 versua 


i Froude 


i Number 


for Bare 








Hull Model 








Ship 








v 


C w/ 


Fr 


Speed 


C w 


C r 


C t 


/c r 






Kts 








% 


% 






xl0 J 


xl0 J 


xlO^ 






0.242 


15.0 


0.204 


0.92 


5.00 


22.2 


4.1 


0.363 


22.5 


0.703 


1.78 


5.53 
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0.363 
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0.329 
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0.218 


0.741 


2.223 
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0.133 


0.744 
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li 
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0.715 


2.473 


5 


ll 


0.255 


0.641 


2.151 


3 


Small 


0.220 


0.549 


2.510 



Construction of Influence Diagrams 

Table IV gives the computed, wave resistance coeffi- 
cients versus Froude number for each of the seven configura- 
tions. Disregarding the small protuberance valves at this 
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time, and dividing each large protuberance coefficient by 
its corresponding bare hull coefficient, Influence diagrams 
showing the effect of the increased sectional area at 
various points along the hull may be constructed for each 
speed. Figure 10 represents the influence diagrams for 
the model, based completely on experimental values. The 
remainder of this section discusses attempts to replicate 
Figure 10 utilizing the linear superposition hypothesis. 

If the techniques applied are successful, the amount of 
model testing and wave record analysis required to obtain 
an influence diagram will be greatly reduced. 

In order to apply the linear superposition hypothesis, 
the PROTUB computer program was modified and adapted to 
the NSRDC 7090 computer. The difference signals were 
shifted one data point at a time with respect to the bare 
hull wave record and the wave resistance coefficient was 
computed for each altered configuration. This is equiva- 
lent to a series of discreet shifts of the protuberance 
along the model hull. Figures 11 and 12 are the non- 
dimenslonalized Influence diagrams derived in this manner 
for the model with large protuberances attached. The 
experimentally-obtained wave resistance coefficients at 
each station are also noted. These curves show that the 
predicted wave resistance coefficients at a shifted position 
are highly dependent upon the Initial position of the 
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FIG 10. 

INFLUENCE DIAGRAMS BASED ON EXPERIMENTAL VSASE 
RESISTANCE COEFFICIENTS 
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STATION SPACING 



INFLUENCE DIAGRAMS EASED ON LINEAR SHIFT OF 

THE DIFFERENCE SIGNALS 
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protuberance and that the predictions do not always agree 
with the actual wave resistance coefficient. 

At this point the physical events which occur when 
the protuberance is moved longitudinally along the model 
during a wave resistance test, which are disregarded in 
the idealized linear superposition theory, should be 
discussed. First, as the probuberance is moved aft along 
the model, the separation between the halves is increased 
and the amount of model interference between the two halves 
is Increased, Secondly, the direction of the streamline 
flow approaching the protuberance differs at each location, 
effectively changing the angle of incidence of the flow. 
Finally, the form of the model varies at each location. 

All of these faotors influence the wave pattern generated 
by the protuberance and could be the causes of the 
differences obtained in the above influence diagrams. To 
minimize these effects, only the protuberance at station 
1 will be used in later shifting applications, since the 
difference signal generated by this configuration most 
closely approximates the protuberance in a free stream. 

Salvesen (5) found that waves tend to die out in 
phase as well as amplitude as they propagate, and conse- 
quently their phase at a distance from the point of genera- 
tion differs from the phase predicted by wave theory. 

Since it is expected that this phenomena is occuring in 
the wave resistance tests conducted for this project, a 
translation was applied to the influence diagram in an 
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attempt to obtain better agreement between the actual and 
predicted wave resistance coefficients. 

By Inspection of Figures 11 and 12 It was decided 
to use a twenty percent aft translation for a Froude number 
of 0.242 and a thirty percent aft translation for a Froude 
number of 0.363. However, for a Froude number of 0.484, 

It was obvious that translation of the Influence diagram 
would be unrealistic, and therefore, no attempt was made 
to apply the technique to this speed. These translations 
yielded wave resistance coefficients at stations 1, 3 and 
5 which are within a few percent of the experimentally 
obtained values. 

To verify this technique, additional model tests 
were conducted with the large protuberance attached at 
stations 2 and 4. At both Froude numbers these tests 
yielded wave resistance coefficients at station 2 which 
were within six percent of the predicted translated 
values, while the wave resistance coefficients obtained 
at station 4 deviated considerably from the predictions. 

The experimentally obtained points have been shown on 
Figure 13. Table V compares the actual and predicted 
resistance coefficients. 
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TABLE V 



Comparison of Actual and Predicted Vave Resistance 

Coefficients 
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at Sta 


Predicted 

C w 
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0.231 


0.218 
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0.362 


0.261 


3S.7 


0.363 
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0.815 


0.741 


10.0 


0.363 
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0.622 


0.715 


13.0 


To 


check the 


validity of 


the linear 


superposition 



technique for short distances, a comparison was made of 
the actual wave resistance coefficients obtained from 
model tests and the wave resistance coefficients predicted 
by shifting the protuberance one station forward or aft, 
with no translation of the Influence diagram. This com- 
parison Is shown as Table VI. It is evident that a higher 
froude numbers a shift of the difference signal of one 
station aft predicts the wave resistance with considerable 
accuracy, while shifting the difference signal forward 
gives no reliable results. 
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TABLE VI 



One Station Shift Comparison for Large Protuberances 
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ft 
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tt 
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2.473 


2.136 


0.607 


24.5 



Data Investigation 

At each Froude number under Investigation, the wave 
record of the protuberance was obtained by the method 
outlined in the introduction. These wave records or 
difference signals are shown as Figures 14, 15, and 16. 

On eaoh figure the assumed starting location, taken as the 
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FIG 16. 

DIFFERENCE SIGNAL^ 
FOR LARGE PROTUBERANCE^ 
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protuberance maximum cross section, of the difference 
signal Is Indicated by^ . 

If the original linear superposition theory were 
completely valid, the difference signals for each protub- 
erance location at a given Froude number would be Identical 
In phase and amplitude, varying only by their starting 
position. This shift In starting position would be directly 
related to the location of the protuberance on the model. 

As can be seen from the difference signal plots, the 
correlation between the curves is not exact. However, 
the following points of agreement should be noted: 

a. All of the difference signals are sinusoidal In 
nature and are quite similar In character to the wave 
records obtained from the model tests. 

b. For a given Froude number, the difference signals 
exhibit fairly good phase and amplitude relationships. 

c. The magnitude of the amplitudes of the signals 
decrease with Increased Froude number, indicating that the 
effect of a protuberance upon the wave resistance decreases 
with Increased speed. 

d. At all Froude numbers the most severe difference 
signal occurs for the protuberance at station 3. 

From this discussion it is evident that direct super- 
position may be possible but that additional analysis of 
the data Is required. 

The amplitude spectra and Fourier sine and cosine 
transforms were plotted so that the data could be compared 
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after It had undergone the first Integration. This process 
gives an Indication of the effect of the deviations between 
difference signals upon the wave resistance coefficient 
and allows comparison of the data after It has been smoothed. 
Figures 17, 18 and 19 show the curves obtained at a Froude 
number of 0.363 for the configurations with large protuber- 
ances attached at stations 1, 3 and 5. The remainder of 
the spectra are Included In Appendix C. 

The agreement among these spectra Is remarkably good , 
particularly at the higher Froude numbers. The resultant 
amplitude spectra and Fourier transform spectra show a 
similar phase relationship at each Froude number as well 
as definite amplitude agreement. 

In addition to the above, plots, the amplitude spectra 
of the difference signal alone were obtained for the large 
protuberance attached at stations 1, 3* and 5 at all three 
Froude numbers. These curves are shown as Figure 20. The 
nature of these spectra Indicates that the difference signal 
obtained Is essentially equivalent to that of the protub- 
erance in a free stream, since In a free stream the protub- 
erance can be represented by a point source, and thus It 
would give a straight line amplitude spectra. 

Although free wave spectra cannot be directly applied 
to shifting methods, their consistent similarity indicates 
that further Investigation of the proposed linear super- 
position technique should be pursued. 
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WAVE COSINE TRANSFORM 

FOR 

LARGE PROTUBERANCE CONFIGURATIONS 
FROUDE NUMBER .303 
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WAVE AMPLITUDE SPECTRA 
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Dlze Comparison 

The effect of protuberance size upon the wave re- 
sistance coefficient was analyzed by comparing the results 
obtained frorajmodel tests with the large and small pro- 
tuberances attached at station 3* The methods of analysis 
were essentially the same methods employed In the 
previous sections and are best explained by the following: 

a. Difference signals, shown as Figures 21 and 22, 
show no dlreot relationship. 

b. Wave amplitude and Fourier sine and cosine trans- 
form spectra, shown for Froude number 0.363 as figures 

23 , 24 and 25 and for Froude numbers 0.242 and 0.434 in 
Appendix C, show relatively good relationships at the 
higher Froude numbers. 

c. Influence diagrams. Figures 26 and 27, are almost 
identical at Fr = 0.434, while there is no agreement be- 
tween the diagrams at the lower Froude numbers. 

d. The ratio of experimental wave resistance coeffi- 
cients at the three speeds are: 

Fr / Cw 

large small 

0.232 0.332 

0.363 1.355 

0.484 1.069 

As is evident from the various curves, the agreement 
between the results is a function of the Froude number 
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with the best agreement occur Ing at the highest speed. 
This is probably because the relative influence of the 
protuberance decreases with Increasing speed. 
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FIG 24. 



WAVE SINE TRANSFORM 

FOR 

LARGE {SMALL PROTUBERANCE CONFIGURATIONS 
FROUDE NUMBER .3G3 
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FIG 25 
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CONCLUSION’S AND K -XO .■PDATIOI-iS 

Contrary to the results obtained by Tggers, Sharma, and 
Ward (1,3) regarding the feasibility of the lateral cut, wave 
slope method of determlng wave resistance, the authors experi- 
enced large variations in the computed wave resistance coeffi- 
cients of individual data record sets. When principles of 
statistical probability were Introduced (Appendix . ) , the 
average wave resistance coefficient used was found to have 
only a 50 % degree of confidence cf being within 6/ of the actual 
value. Therefore, all results in this paper must be viewed 
with a certain degree of uncertainty. 

At the present time, influence diagrams can c-e obtained 
only by extensive model testing. Attempts to accurately 
produce an Influence diagram by applying the pinciple of linear 
superposition of wave records were met with varying degrees of 
success. The following positive trends were observed: (a 1th 
increasing broude number the relative effect of the protuber- 
ance upon the wave resistance decreased, yielding better 
correlation between the experimental and predicted wave resist- 
ance coefficients. (b) The linear superposition technique of 
developing influence diagrams appears to be valid only over 
small Intervals. (c) At the lower Troude number, translation 
of the influence diagram was required to obtain any decree of 
correlation between experimental and predicted wave resistances. 
This technique produced relatively good agreement within an 
interval of two stations on the model. 
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Then two protuberances of different size were tested at 
the same location, no discernable agreement or uniform trends 
were found at the lower p’roude numbers. Only at the highest 
"roude number, where the effects of the protuberance are 
diminished, was any consistency observed. 

Based on the foregoing results, the following recommenda- 
tions are made for future work In this field: (a) The statis- 

tical nature of the data should be investigated to determine 
the approximate number of data records required to obtain 
sufficient confidence in the final wave resistance values. 

Phis appears to be especially true for the type of wave slope 
measurement used in the present investigation. (b) oince 
there have been no previous wave resistance tests on Festroyer- 
type models, the percentage relationship between the wave and 
residuary resistances is not known. It is possible that the 
low values obtained in the present tests reflect an error in the 
method or instrumentation, Further equipment checks in both 
areas should be made. (c) Curve fitting techniques should be 

r 

applied to the experimentally -obtained wave resistance coeffi- 
cients to improve the quality of the influence diagrams. (d) 
Theoretical investigation into non-linear wave theory should v e 
conducted to discover, if possible, simple adjustments to the 
linear superposition principle so that acceptable agreement 
between the experimental and predicted wave resistances and 
between results for different-sized protuberances can be obtained. 
The similarities between the difference signals and the wave 
spectra noted earlier encourage this effort. Success would 
allow one to construct a valid influence diagram for wave resist- 
ance which does not seem possible using the present linearized 
approach . 
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IX A 



’io: 



T^ 



;lo; 






A successful wave slope measurin' probe based on the 
resistance wire principle was developed in 1965 for use at 
Tebb in conjunction with the 2KC excitation from a two-channel 
brush .j -4622 -00 amplifier and Mark II recorder. The 
probe and the calibration scheme is shown in figure A-l, the 
corresponding circuit diagram is shown in "igure A -2, and 
photographs of the equipment in use are shown as figure A -3. 
'"’he scheme is an adaptation of the wave height measuring pro^e 
previously developed in connection with this effort (5). 
and employs three parallel stainless -steel wires of .004 
inch diamete~, spaced at 3/4 inch intervals between care- 
fully drilled holes in 1 lexiglass blocks held by a T-shaped 
frame. "he frame is mounted on a slider 'which allows known 
variations of vertical height of + 3" at approximately 
y" intervals and this in turn is on a plate which can be 
rotated in either direction to produce a slope variation of 
.20 radians In steps of 0.10 for the slope calibration. 

"he height variations are used to calibrate wave height 
probes mounted on the same frame or in the present case to 
check the insensitivity of the wave slope probe to changes 



in wave elevation. 



4*. ""his appendix has been extracted from " 'xperimental 

. eterraination of 'ave ’ esi stance of a chip odel from 
lateral ave Hope Measurements" by I. r . hard, dep. 19-7 
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As shown In Figure A -2, the measuring circuit consists 
of a Wheatstone bridge with two of the resistances and 
R^ replaced by the equivalent path resistance of the water 
between each set of wires AB and AC. This scheme Is known 
In the case of the height probe to give a good linear response 
to changes In wave height over a range of several Inches 
when the ratio of resistance: 

r 3/R4 = R 2/R]j = K (A-l) 

Is properly chosen; In this case the value of K Is approxi- 
mately 20 based on the equivalent D. C. resistance of the 
water elements. In the case of the slope probe the elements 
are used differentially, as shown. This circuit Is more 
symmetric and actually easier to balance than the height- 
probe circuit where an electrloal resistor must be used to 
attempt to reproduce the water effect. A combined resistance 
and capacitance balancing circuit Is provided within the 
bridge circuit to assist that available In the amplifier. 

As mentioned previously, good linear response has been 
established for this probe in operation over the range of 
t 0.20 radians ample for the usual model wave system and 
relatively good insensitivity with changes In wave height 
over the range of 1 l|r Inches required for the present tests. 
Care is taken to keep the wires clean and at constant tension. 
No difficulty with this arrangement has been noted due to 
variations In surface film, nor does the exact spacing of 
the wires seem to be a critical factor as has been noted by 
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others. This might be a fortuitous result of the relatively 
high exciting frequency of 2KC provided by the Brush amplifier. 

The foregoing statements refer to the characteristics 
noted under static conditions and do not include any effects 
due to dynamic action, such as meniscus oscillation at the 
frequency of the waves being measured. Although Pearlman 
ran a series of such tests on a similar device and found 
no Important effects at the frequencies of importance, the 
present device should be subjected to a similar series of 
tests for the sake of completeness. 
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APPENDIX B: COMPUTER PROGRAM 



A basic program computerizing Equation (14) Is given 
by Ward (3). An extensively modified version of It was 
adapted for this paper, although the fundamental computa- 
tions have remained Intact. In addition to computing a 
wave resistance coefficient, the program (named "PR0TU3") 
developed for this project can subtract a bare hull record, 
point by point, from a wlth-protuberance configuration 
record, and compute the difference signal wave resistance 
coefficient. Additionally, It can shift this difference 
record forward or aft along the model, multiply each slope 
value by any desired constant, add the modified signal back 
to the bare hull record, and determine the wave resistance 
coefficient of the altered data. 

PROTUB, as listed In Figure 3-1, is written In GE 
time-sharing fortran. The input data consists of NZETA J 
wave slope values which are stored and referenced In the 
SFILE statement. Eaoh data point Is divided by the calibra- 
tion faotor for that configuration (CALI or CAL2) to obtain 
the ZETAY(N) set in radians, which are the proper units 
for processing. 

Equation (14) may be re-wrltten In series form as: 



5. Capitalized words correspond to names as programmed In 




(3-1) 
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PROTUB 



The Incremental value of L{ Is given In Ref (3) as 

^ = ^/Jobey CB-2) 

where h Q ff is a hypothetical effective tank width. If 

it Is assumed that k Q b ef » f = loTT , then ,4 ^4 =— — = 0.10 

NIETA >0rr 

3 y = SY =^1 ZETAY(h) S!H(vjA 0 x) AA 

A/- / 

V i£7^ _ . .S . . , . ... { (3-3) 



Oy = C* = 



ZETAY(h) CO s(wA„x) AX 



V s / 



6. This simplifying assumption may be justified by consider- 
ation of the tank geometry. The actual effective tank 
width considering the model offset equals 11.38 ft. 



b ef f 



rr 





Ft 



Therefore b e ff is generally equal to or greater than 
the actual value, giving a about equal to or smaller 
than required from Equation (3-2) . This provides a 
finer integration than is required by Equation (3-2). 
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where : 



X = DELX = 7i'i( model speed In ft/sec) " PAPE?, 
(interval between readings in mm) / VP (paper 
speed in mm/sec) 

Therefore DELX equals the distance between data 
points in feet 



X = N * DELX, N = 1,2,3 NZSTA 

w is a function of J and is computed from Equa- 
tion (11) for each Increment of U . 

Letting CAYO = k Q and factoring DELX out of the expression, 
Ni£TA 

5 u = C -3ETAY (N) * SiN C W * CAYO* Kl * D£o0 

* V-I 



Cy - ^ A Z£TAYW)*COS(W* CAY0*N* DEoOj 



(3-4) 



Equations (3-4-) solve the Fourier transforms. How- 
ever, the trigonometric identities, 

sin(X + A X) « sinX cosAX + cosX slnAX 
cos(X *A X) = cosX cosAX - sinX slnAX 
are utilized in P30TU3 to reduce computer core time. 
Letting kobgf ~ C.mY03 — iotr 



S = WETSR = model stir face -ft 

o ^ 

°w = 2 ( DSLX) ^ 

WETSR- *CA YOB 

j=o. 



el su 
£ 



SY 2 + CY 2 

w^*(2w 2 «l) 



(3-5) 

where the entire expression's doubled because only | the 
wave pattern is recorded by the probe. Note that the final 
resistance coefficient is the product of a constant and the 
wave amplitude spectra, squared, times a weighting function, 

_1 J , which decreases quite rapidly with Increase in 

w^(2w 2 -1) J u or w. 
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APPENDIX C: WAVE SPECTRA PLOTS 



Tha remaining amplitude and Fourier sine and cosine 
spectra, not Included In the body of the report, are shown 
In this section. The curves were plotted by the NSRDC 7090 
computer. Due to the limitations of the curve plot sub- 
routines used, the following notation was adopted for identi- 
fication of the curves: 

First letter denotes: 

A Bare hull configuration 
B Large protuberance at station 1 
C Large protuberance at station 3 
D Large protuberance at station 5 
E Small protuberance at station 3 

First number denotes: 

1 Froude number equals 0.242 

2 Froude number equals 0.363 

3 Froude number equals 0.484 

AVE denotes that the curves are plotted from data averaged 
over three model tests. 
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APPENDIX D: ERROR ANALYSIS 



Data Inconsistency Record 

Although the three data tapes retained for each configu- 
ration and speed had an average speed within t 0.01 ft/seo, 
the wave slope records were not Identical, and significant 
phase shifts developed toward the middle of most of the 
records. To investigate this source of error four sets 
of data were artibrarlly selected. These were: 

Set Speed Configuration 



242 


Large protuberance 


at 


station 1 


363 


Bare Hull 






363 


Large Protuberance 


at 


station 2 


484 


Bare Hull 







When individual tapes of each set were analyzed, the 
extreme wave resistance coefficient had an average devia- 
tion of almost 10 % from the mean. The deviation amoung the 
four sets of data so analyzed are shown in Table D-I. 

TABLE D-I DATA DEVIATION 







o 

o 


Individual Data Tapes 


Set 


Tape 1 


Tape 2 


Tape 3 


c w 


Max Dev from C, r 
w 


1 


.3129 


.3545 


.3114 


.3290 


7.8 


2 


.6 755 


. 7780 


. 6816 


.7033 


10.7 


3 


.8404 


.7394 


.6 790 


.7529 


11.6 


4 


2,348 


2.382 


2.047 


2.244 


8.8 



Average 9 . 7 



101 * 



There Is strong reason to believe that the wave records 
are statistical In nature, and that many more data tapes 
are required to attain a suitable degree of confidence in 
the average value. Applying the student ’T* distribution 

- 2 

for small samples, If C w and s are the mean and variance of 
a sample of size N from an assumed normal distribution, 
nC^.<^“) . where^ and <T are unknown, then the percent 
deviation of C w from the mean with a given confidence Interval 
may be readily determined. 

Each of the four sets of data Individually investigated 
and compared in Table D-I experienced considerable variation. 
Therefore, probable deviations of from the theoretical 
mean were determined for each set at four confidence levels 
and averaged. These steps are tabulated in Figure D-l, 
which shows the percent deviation of the accepted C from 

w 

the mean versus its respective degree of confidence. For 
comparative purposes, the same values of s and are used 
to calculate the confidence levels were ten samples to be 
taken. It is quite apparent that, from a statistical view- 
point, three runs per set do not produce sufficiently reliable 
values of wave resistance coefficients. 

Calibration Factor 

The Brush recorder and amplifier exhibited some fluctua- 
tion and non-linearity during testing. A sample calibration 
factor plot for the bare hull configuration is shown in 
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DEVIATION OF C w FROM MEAN IN PERCENT 



FIGURE D-l CONFIDENCE LEVELS FOR WAVE RESISTANCE COEFFICIENTS 



SET 


c w S/JJ 


% DEVIATION FROM 


MEAN AT <jo CONFIDENCE 
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1 


.3290 .0143 


3.55 


8.2 


12.7 


18.7 


2 


.7033 .0337 


3.92 


9.0 


14.3 


20.7 


3 


.7529 .0471 


5.20 


11.8 


18.3 


27.0 


4 


2.244 .0830 


3.02 


7.0 


10.8 


15.9 


AVERAGE DEVIATION 


3.90 


9.0 


14.2 


20.5 




DEGREE OF CONFIDENCE IN PERCENT 
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Figure D-2. Depending on the manner of securing the probe 
after each rotation when calibrating, an error of - i mm. 
at .2 radians was observed. This gives an expected calibra- 
tion error of t 1.25 mm/rad. 

As the probe Is raised or lowered with no rotation, 
a perceptible soale deflection was detected. However, the 
resulting effect on the wave record tends to be self -com- 
pensating and this factor Is mentioned for Information only. 
Therefore, since the average calibration factor Is 

about 95 mm/rad and Is squared In the calculations, per- 

2 

centage error Is approximately (± 1.5/95) ^ 2|$. 

Wave Record Reading Error 

Because the human element Is Involved In reading the 
data tapes, a reading error of measurable magnitude Is likely 
to be Incurred. To Investigate this problem, each of the 
tapes In the fourth set of data mentioned above was read 
individually by each of the, authors and separate values of 
computed. Table D-H summarizes the results. From this, 
a gross reading error of * 1 3/k% Is anticipated. 

Table D-II Reading Error 

Protuberance located at station 2, Froude number = 0.363 



Reader 


1st Tape 


2nd Tape 


3rd Tape 


Average 


FJM 


.8337 


.7216 


.6732 


.7278 


CAB 


.8383 


.7648 


.6712 


.7500 


AJJ 


.8492 


.7317 


. 6926 


. 7448 



Average for all readings = .7409 



Maximum deviation from average = 7409-7278 

.7409 = 



Figure d-2 bare hull calibration factor 
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Overall Sensitivity of Method. 



If a sufficiently large sample of wave records were 
taken, the percent deviation of C w at 90 $ confidence could 
readily be reduced to a few percent. Maximum combined 
calibration and reading errors constitute about 2^+1 3/4 
or 4$ error. Therefore, the longitudinal cut method, when 
conscientiously followed, with an adequate number of records, 
should provide wave resistance values accurate to within 
6$ or less. However, when only three samples are used, a 
6 % accuracy is available with only a 50$ confidence level. 
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